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Abstract 



> 

^^ ■ For the complex Ginzburg-Landau equation on a large periodic interval, 

(<— ^ ■ we show that the transition from defect- to phase-turbulence is more accu- 

On . rately described as a smooth crossover rather than as a sharp continuous 

■^: 

J>^. transition. We obtain this conclusion by using a powerful parallel computer 

^ . 
Q . to calculate various order parameters, especially the density of space-time 

,^ . defects, the Lyapunov dimension density, and the correlation lengths of the 

field phase and amplitude. Remarkably, the correlation length of the field 

amplitude is, within a constant factor, equal to the length scale defined by 

the dimension density. This suggests that a correlation measurement may 

suffice to estimate the fractal dimension of some large homogeneous chaotic 

systems. 
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Recent advances in laboratory technique 0] and in computer simulation 0-§[ have 
opened up the study of boundary- independent spatiotemporal chaos in large homogeneous 
sustained nonequilibrium systems 0. Many fundamental questions remain unanswered 
about such chaotic systems, e.g., what different states can occur, how transport depends on 
different states, and what kinds of bifurcations separate one state from another. An espe- 
cially interesting question is whether ideas from statistical mechanics might be applicable 
to chaotic nonequilibrium systems in the thermodynamic limit of infinite system size [0-^ . 

A significant step towards understanding some of these questions was recently reported 
by Shraiman et al [§]. These researchers studied different spatiotemporal chaotic states of 
the one-dimensional complex Ginzburg-Landau equation 

dtu{x,t) =u + {l + ici)dlu-{l-ic3)\u\\, xg[0,L], (1) 

on a large periodic interval of length L = 1024, which they assumed to be large enough to 
approximate the thermodynamic limit of infinite system size. Here the variables t and x 
denote time and position respectively, the complex-valued field u{x,t) = pe*"^ has magni- 
tude p{x, t) and phase 0(x, t), and the parameters Ci > and C3 > are real- valued. Eq. (|Tp 
is an important model of spatiotemporal chaos because it is simple, experimentally relevant 
T0| and universal P: spatially extended systems that undergo a supercritical Hopf bifur- 



cation from a static to oscillatory homogeneous state reduce to Eq. ([l|) sufficiently close to 
the onset of the bifurcation. Interesting dynamics are predicted and are observed beyond 
the Newell line C1C3 = 1 since all plane wave solutions of Eq. (||) are linearly unstable to the 
Benjamin- Feir instability for C1C3 > 1 0. 

Shraiman et al summarized their simulations in the form of a phase diagram in the 
C1-C3 parameter plane (Fig. 3 of Ref. [Q]). Based mainly on calculations of the density of 
space-time defects n£, |Tl|, this diagram showed continuous and discontinuous transition 



lines (analogous to second- and first-order phase transitions) separating defect-turbulent 
from phase-turbulent states |TT|. Of special interest to us is the continuous chaos-to-chaos 



transition line Li in their Fig. 3, which occurs for Ci > 1.8. It is somewhat mysterious why 



the density n£) decreases to zero at an Li line that is distinct from the Newell line: in the limit 
of infinite system size and of infinite time, what prevents defects from forming anywhere to 
the right of the Newell line (ca > 1/ci)? The mystery of the Li line can be partly appreciated 
by trying to reason by analogy to equilibrium statistical physics. Assuming that the chaotic 
fluctuations of Eq. (|I|) act as a finite-temperature ergodic noise bath and observing that 
the derivatives in Eq. (|l]) represent short-ranged interactions between different parts of the 
field M, we would not expect distinct phases at finite temperature in one-space dimension 
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Because so little is known about possible critical phenomena of large homogeneous 
nonequilibrium systems and because Eq. (|I]) is such an important model, we have tried 
to characterize more carefully the dynamics near the Li line for the fixed parameter 
value Ci = 3.5. By calculating various order parameters over length scales as large as 
L < 10^ and over time scales as large as T < 10"^, we are able to show below that the 
change from defect- to phase-turbulence near the Li line is more accurately described as 
a smooth crossover rather than as a sharp continuous transition with power-law scaling of 
order parameters 0. It is then a possibility that phase turbulence (i.e., a chaotic state 
with n/) = in the thermodynamic limit) does not exist although we can not settle this 
with our present computer resources. Our calculations also confirm several points predicted 



by Shraiman et al and by other researchers |13|, e.g., that the spatial correlation function 



of the phase should decay exponentially inside the phase-turbulent regime. We have also 
studied whether the dimension density 5 (Lyapunov fractal dimension per unit volume) is a 
useful order parameter for characterizing changes in spatiotemporal chaotic states P,|3|,p!4|. 
The dimension density defines a dimension correlation length ^s = S^^^^ p which is the 
characteristic size of dynamically independent subsystems of spatial dimensionality rf 0. A 
comparison of C,5 with other characteristic length scales as a function of the parameter c^ 
gives the remarkable result that C,s is, up to a constant factor, equal to the spatial correlation 
length of the field magnitude, C,p, from the Newell line to beyond the Li line. 

An important resource for the advances reported below was a CM-5 parallel computer 



T5| , which facihtated the study of much larger space-time and parameter regions than 
were previously conveniently accessible. Our numerical methods for integrating Eq. (|I]) 
and for calculating related order parameters are similar to previous methods except for 
modifications of algorithms and codes to take advantage of the CM-5's scalable parallel 
architecture ||^. For most of our simulations, we used a time step At = 0.05 and a spatial 
resolution of two Fourier modes per unit length; this space-time resolution was dictated 
by the need to detect isolated space-time defects when estimating the density rir, |^. Our 
initial condition for most runs was u{x^ t = 0) = 0.4 + 7] where rj was a uniformly distributed 
(5-function correlated random variable in the interval [—0.02,0.02]. We typically averaged 
over 64 or more random initial conditions spanning an integration time of 2 x 10^ so that 
the total effective integration time was perhaps as long as 64 x (2 x 10^) = 1.3 x 10^ time 
units. 

Before discussing our results, we note that for ci = 3.5, for an integration time of 
T = 10^, and for a periodic interval of length L = 1024, Shraiman et al argued the existence 
of the Li line using two key observations 0: (1), that the density no vanished as a power 
law n£, oc (c3 — C3)" with exponent a ^ 2 and with C3 ~ 0.77 > c^^^^'-'- = \jc\ = 0.286; 
and (2), that the correlation time r of phase fluctuations ITO diverged as a power law 



also at C3, as the inverse of the defect density, r oc l/nc- For C3 < C3, Shraiman et al 
observed a less-disordered phase-turbulent regime with ud empirically equal to zero and 
with slower-than-exponential decay of temporal correlations 0. If defects do not occur in 
the thermodynamic limit, a perturbation theory in the small quantity e = C1C3 — 1 yields a 
simpler description of phase turbulence near the Newell line, e — > 0. In that limit, Eq. (|l]) 
reduces to the Kuramoto-Sivashinsky (KS) equation [|1^, 



d,(P = -edl<P - ]^c\{l + cl) d^J - (ci + C3) (a,0)2, e = C1C3 - 1, (2) 

and the amplitude p becomes an algebraic function of a spatial derivative of the phase, 

p ^ 1 - |9^0. (3) 



Some of our calculations below provide the first quantitative comparisons of phase turbulence 
as described by Eqs. (0) and (|^) with phase turbulence as empirically observed in Eq. (|1]). 
For the parameter value ci = 3.5, for a system size L = 4096, and for an effective 
integration time of T = 10^ (after allowing transients of duration 10^ to decay), we find in 
Fig. |I| that no is finite substantially to the left of the Li line as calculated in Ref. 0. Far to 
the right side of the Li line, our data in Fig. 0(a) approximately reproduce the previously 
reported P] power-law scaling with exponent a ^ 2. Closer to the Li line, a least-squares fit 
of the three-parameter expression a(c3 — C3')" to the nine left-most points gives a much larger 
exponent a ~ 6.8, with an onset of phase turbulence {ud = 0) at 63' = 0.74 < C3 = 0.77. 
Assuming equal errors bars on all data points, we find the chi-square value for the fit to be 
X^ = 4.6 X 10^^^. The increase in the exponent with increased space-time resolution suggests 
that a power-law scaling is inappropriate. As shown in Fig. |l](b), we find a better fit of the 
same data with the functional form 

nD = a exp (-6/ (cg - C3")") > (4) 

which is the expected behavior for thermodynamic Gaussian fiuctuations of the phase gra- 
dient dx(f) if large values of the latter are the reason for defect nucleation 0. If we set a = 1, 
a least-squares fit of Eq. (^) to the nine left-most data points yields the three parameter 
values a = 0.66, b = 0.98, and C3" = 0.70 < C3' = 0.74 with a x^ = 8.1 x lO^^^^ jf ^g gg^ 
C3" = c^^"^^'' = 0.286 to test whether Eq. (^ is consistent with the onset of phase turbulence 
at the Newell line, a least-squares fit (again to the nine left-most points) gives the parameter 
values a = 0.018, b = 0.017, and a = 8.8 with a substantially poorer x^ = 8-2 x 10^^^. Our 
data spanning the crossover region evidently lie too far to the right of the Newell line to 
determine whether the defect density goes to zero before or at this line. 

To test independently the important implication of Fig. |1| that a crossover occurs (so that 
phase turbulence may not be a distinct phase from defect turbulence), we have calculated 
other order parameters over the same parameter range. The occurrence of a crossover is 
supported by Fig. |, which summarizes correlation times r of the phase (f){x, t) |]l6l as a 



function of C3. Fig. ^(a) shows that r does not diverge to infinity at the Li hne (as reported 
in Ref. 0) but instead is large and finite a httle bit to the left of the Li line. We find 
that time correlation functions of the phase decay exponentially down to at least C3 = 
0.77 but do not have a simple functional form for smaller C3 values. Recent arguments 
and calculations have been made |T^ that time correlation functions should decrease as a 



stretched exponential for the phase-turbulent regime of Eq. ([l|) in which case a correlation 
time can not be meaningfully defined. However, the change from exponential to stretched 
exponential behavior with decreasing C3 has not yet been carefully studied. Shraiman et 
al conjectured that the critical scaling of r at the Li line followed from the scaling of ud 
with T oc l/riD. Fig. 0(b) shows that this is approximately correct to the right of the Li line 
but breaks down closer to this line since the product ttij:, is no longer constant: the defect 
density goes to zero faster than the correlation time increases. 

In Fig. 0(a), we have calculated the phase spatial correlation length $,^ |16| as a function 



of C3. Shraiman et al argued that C,(f, should be finite in the phase turbulent regime of Eq. (|I]) 
and estimated its value indirectly by calculating a phase diffusion coefficient D = 1/^^ from 
phase-gradient correlations [^ . Exponential decay of spatial correlations is also expected for 
phase turbulence if the latter is described at long-wavelengths by the Kardar-Parisi- Zhang 
(KPZ) Langevin equation [|l^. By going to quite large system sizes {L = 10^) and to long 



integration times, we have verified directly that the phase spatial correlation function ||T^ 
decays exponentially well to the left of the Li-line as shown in Fig. 0(a). As the parameter C3 
decreases, the quantity ^^ varies smoothly through a local maximum near the Li line, and 
then increases steadily until we can no longer estimate its value accurately with our computer 
resources. The smooth variation of ^^ through the Li region is consistent with a crossover 
rather than with a sharp transition. The apparent divergence of ^^ upon approaching the 
Newell line, e — > 0, can be understood semiquantitatively as shown in Fig. |^(a) by a scaling 
argument [|I^ that predicts ^^ oc e~^/^. The agreement is within about 10%. 

The phase correlation length C,^ is the same as that of the field u itself [Q , but there is 
a separate, generally shorter, correlation length scale C,p associated with fiuctuations of the 



field amplitude p (also with the phase gradient dxcf)). Fig. H(b) compares the reciprocals of 
the phase and amplitude correlation lengths with the Lyapunov dimension density 5, whose 
reciprocal defines the dimension correlation length ^s discussed above [Q. Up to constant 
factor of 1.4, we find that the amplitude correlation length equals the dimension correlation 
length ^5 over a substantial range of parameter C3. (An independent and related result was 
also recently reported by other researchers [|19|.) This remarkable result suggests that the big 



fractal dimension of some large homogeneous chaotic systems might be accurately estimated 



we 



by simple correlation function calculations. In work that we will discuss elsewhere |I4 
have also calculated the variation of C,5 across a nonequilibrium Ising transition in a two- 
dimensional coupled map lattice p. Although the agreement is not quite so striking, ^s still 
matches closely the correlation length associated with fluctuations of the magnitude of the 
site variables. 

In Fig. ^ we make two final comparisons of how phase-turbulence, as described by 
the adiabatic approximation Eq. (^ and by solutions of the KS-equation Eq. (Q), agrees 
with numerical solutions of Eq. (|T]). The dimension density 6 of the KS-equation has been 
calculated to be 5 = 0.230 for the rescaled parameterless version of the KS-equation [ED 



(9i0 = —9^0 — d^cp — 09a;0. Restoring the original space, time, and magnitude scalings gives 
the following ci and C3 dependence of the dimension density for KS phase turbulence: 



'2(ciC3 - 1: 



1/2 



In Fig. §(a), we compare Eq. (^) with our empirically determined values of 6 for Eq. (|l]) 
from Fig. ^(b). The agreement is good up to about C3 = 0.5 (e = .75) and then there is an 
increasing deviation of the actual solutions from Eq. (^. This deviation with increasing C3 
may arise because the adiabatic approximation Eq. (^ breaks down or because higher-order 
terms in the KS-equation are renormalizing the dimension density. Fig. ^(b) gives some 
further insight by comparing the mean-square fluctuation of p from Eq. (|1|) with the mean- 
square fluctuation of p as given by Eq. (^). We observe a previously unreported power-law 
scaling of these amplitude fluctuations with exponent a = 4 from the Newell line to near 



the Li line. Sufficiently close to the Newell line, an exponent of 4 is predicted by rescaling 
the solutions of Eq. (Q). The adiabatic approximation is satisfied over a larger range in C3 
than the agreement between dimension densities. 

In conclusion, we have used a powerful parallel computer to characterize more carefully 
the change from defect- to phase-turbulence near the Li line in the periodic one-dimensional 
Ginzburg-Landau equation in the limit of large system size. Instead of a sharp continuous 
transition with power-law scaling of order parameters [@], we found a smooth crossover with 
new and anomalous structure near the Li line, e.g., the variation of ^,4, in Fig. |^(a) and the 
change in slopes of dimension density and amplitude fluctuations in Fig. ^. We confirmed 
recent predictions P,p!3| that the phase correlation function decayed exponentially well to the 



left of the Li line, with the related correlation length being finite and large. We also found 
two length scales associated with the field u, a long scale associated with phase fluctuations 
and a short scale C,p associated with amplitude fluctuations. Surprisingly, the length ^p 
equals, up to a constant factor, the dimension correlation length ^s associated with the 
dimension density. This suggests that spatial correlations of certain observables may suffice 
to estimate big fractal dimensions of some large homogeneous chaotic systems. 

Numerous interesting questions remain for future study. Our calculations leave open the 
theoretical question of whether phase turbulence {no = 0) exists in the thermodynamic limit. 
More generally, it is still not known whether a chaos-to-chaos nonequilibrium transition 
can occur in an infinite one-dimensional system. Some of our results might be tested by 
experiment [|10|, e.g., the variation of phase correlation length (Fig. 0(a)) and the scaling of 



amplitude fluctuations near the Newell line (Fig. |(b)). It would be interesting to extend our 
calculations to other parts of the C1-C3 parameter plane, e.g., to understand the hysteretic 
bichaotic regime Ci < 1.8, for which defect-turbulent and phase-turbulent states evidently 
coexist in parameter space 0,|]. If phase-turbulence is just a small ud limit of defect- 
turbulence, it is more difficult to understand the coexistence of different states of identical 
symmetry. Finally, it would be interesting to repeat similar calculations in two- and three- 
space dimensions, for which the point-like space-time defects in one-space dimension are 
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replaced by long-lived topological defects such as vortices and lines [^ . 
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FIGURES 
FIG. 1. (a) Log- log plot of the space-time defect density no versus the distance C3 — C3' to the 
fitted point C3' where the density goes to zero (onset of phase turbulence) for system size L = 4096, 
integration time T = 10^, and an average over 64 randomly specified initial conditions. The arrow 
labeled "Li" indicates the position of the Li-line for parameter value ci = 3.5 0. The smallest np 
value corresponds to a count of 200 defects. The two solid lines were drawn to indicate the previous 
and present best estimates of the exponent q of a power-law scaling. The crosses are the data from 
Ref. P]. In (b), we find that Eq. H with exponent a = 1 gives a better fit of the same data, with 
an onset of phase turbulence at C3" = 0.70. The straight line is a plot of Eq. HI over the range of 
its fit. 

FIG. 2. (a) Plot of the reciprocal square root of the correlation time r (as estimated from 
the asymptotic exponential decay of time correlation functions of the phase e**^ = m/IuJ) versus the 
parameter C3. The arrow indicates where the Li-line occurs for ci = 3.5. We used a system size 
of L = 4096, an integration time of T = 10^, and an average over 64 random initial conditions. 
The crosses denote the similar data from Ref. Q. (b) For the same numerical parameters, a plot 
of the product of correlation time with defect density, tud, showing that the scaling r ex I/hd 
breaks down near the Li-line. 

FIG. 3. (a) Plot of the phase correlation length ^^ for solutions of Eq. || for system sizes of 
up to L = 10^, integration times of up to T = 2 x 10^, and averages 64 randomly chosen initial 
conditions. The crosses denote the similar data from Ref. Q. The solid curve is the analytical 
expression obtained by scaling the finite correlation length of the parameterless KS-equation ||lH| . 
(b) Plot of the Lyapunov dimension density S P and the reciprocals 1/^p and 1/^^ of the amplitude 
and phase correlation lengths. The reciprocal length l/£,p (open circles) has been scaled by a 
constant factor of 0.7 to emphasize the close agreement with 5. The positions of the Newell- and 
Li-lines for ci = 3.5 are denoted by the arrows labeled "N" and "Li" respectively. 
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FIG. 4. (a) Comparison of the dimension density 6 for solutions of Eq. (|T]) with the rescaled 
dimension density of the KS-equation, Eq. (^), for ci = 3.5. (b) Comparison of the mean-square 
fluctuations of the amplitude p as calculated from the Id CGL equation and as calculated from the 
adiabatic approximation, Eq. (|3|), with cp determined from Eq. (|l]. In both (a) and (b)), the arrows 
labeled "N" and "Li" denote the positions of the Newell- and Li-lines respectively for ci = 3.5. 
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